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Very  fast  procedures  for  generating  normal  variables  may  be  based 
on  representing  the  density  function  as  a  mixture,  with  the  dominant 
terms  chosen  so  that  they  lead  to  fast  procedures  in  the  computer.  See 
[l]  -  [3] •  There  is  always  the  problem  of  handling  the  tall  of  the 
normal  distribution.  Variables  from  the  tail  are  needed  so  infrequently 
that  convenience  is  a  more  important  consideration  than  speed  in 
searching  for  methods.  The  following  method  is  very  convenient  -  easy 
to  understand  and  easy  to  program.  At  the  same  time,  it  is  reasonably 
fast,  requiring  a  logarithm  and  a  square  root. 

The  idea  is  to  transform  the  tail  of  the  normal  distribution  to 
the  unit  interval  and  then  use  the  rejection  technique.  It  goes  as 
follows  s 


To  generate  a  standard  normal  variable  X,  conditioned  by  X  >  a, 
-5x* 

i.e. ,  with  density  ce  ,  x  >  a,  generate  pairs  of  uniform  over 
(0,1)  random  variables  U-^Ug  until 


(1)  D2  <  a(aZ-2  In 

then  put  X  =  (a2-2  In  U^. 


The  density  of  U^,  given  condition  (l),  is  a  multiple  of 

(a2-2  In  unr*  0  <  u..  <1,  hence  the  density  of  X  =  (a2-2  In  U.)^ 

-.5x2 

is  a  multiple  of  e  ,  a  <  x  <  a> . 
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Graphs  of  Ug  =  a(a^-2  In  u^)  for  a  =  1  and  3  are  plotted 
in  this  figure.  When  a  =  3,  the  probability  of  event  (1)  is  .88, 
so  that  the  efficiency  of  the  rejection  technique  is  satisfactorily 
high. 
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